Martingale Models for Quantum State Reduction 

S. L. Adler* 1 , D. C. Brody t2 , T. A. Brun* 3 and L. P. Hughston* 4 

"Institute for Advanced Study, Einstein Drive, Princeton, NJ 08540, USA 
^Blackett Laboratory, Imperial College, London SW7 2BZ, UK 

f^) . ^Department of Mathematics, King's College London, The Strand, London 

<N ' WC2R 2LS, UK 



o 



Abstract. Stochastic models for quantum state reduction give rise to statistical laws 
f^) | that are in most respects in agreement with those of quantum measurement theory. 

Here we examine the correspondence of the two theories in detail, making a systematic 

use of the methods of martingale theory. An analysis is carried out to determine 

£> ' the magnitude of the fluctuations experienced by the expectation of the observable 

C*~) ■ during the course of the reduction process and an upper bound is established for 

the ensemble average of the greatest fluctuations incurred. We consider the general 
r — ' projection postulate of Liiders applicable in the case of a possibly degenerate eigenvalue 

^^ . spectrum, and derive this result rigorously from the underlying stochastic dynamics 

for state reduction in the case of both a pure and a mixed initial state. We also analyse 
the associated Lindblad equation for the evolution of the density matrix, and obtain 
an exact time-dependent solution for the state reduction that explicitly exhibits the 
transition from a general initial density matrix to the Liiders density matrix. Finally, 
we apply Girsanov's theorem to derive a set of simple formulae for the dynamics of 
the state in terms of a family of geometric Brownian motions, thereby constructing an 
rz_| explicit unravelling of the Lindblad equation. 
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1. Introduction 

According to nonrelativistic quantum mechanics, the evolution of the state of an isolated 
quantum system is described by a deterministic unitary transformation, governed by the 
Schrodinger equation. The behaviour of the state of a quantum system following the 
measurement of an observable is less well understood, however, and has been the subject 
of much debate. In quantum measurement theory it is usually assumed that when the 
measurement of an observable with a discrete spectrum is carried out on a system 
prepared in a prescribed initial state, then the state reduces randomly to one of the 
eigenstates of the observable being measured. This is the so-called projection postulate 
of von Neumann []I[ , which was later generalised by Liiders |2| to handle the case of a 
measurement of an observable with a degenerate spectrum. The Liiders postulate has 
the virtue of being unambiguously applicable whether or not the initial state is pure, 
and whether or not the eigenvalue spectrum is nondegenerate. 
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Stochastic extensions of the Schrodinger equation have been increasingly attracting 
attention as plausible dynamical models for state vector reduction in quantum mechanics 
H |], ||, H [7], H, 3 . In such models, the Schrodinger equation is generalised to take the 
form of a special type of stochastic differential equation on Hilbert space. For example, 
given a Hamiltonian H and a commuting observable F with a discrete spectrum, there 
exists a natural stochastic differential equation generalising the Schrodinger equation 
with the property that starting from a given initial pure state, the system evolves 
randomly in such a way that asymptotically it reaches one of the eigenstates of F 
with the correct quantum probability. More generally, given a compatible family of 
observables F a (a — 1, 2, • - • , r), each of which commutes with H, then a similar result 
holds, with an asymptotic reduction to one of the common eigenstates of the given 
family of observables. For recent reviews outlining the development of this approach, 
with extensive references, see (JTOj] and ||T] . 



The purpose of this paper is to analyse in some detail the statistical laws associated 
with stochastic extensions of the Schrodinger equation, and to show in particular how the 
projection postulate in the general form due to Liiders can be derived as a consequence 
of the dynamics. The plan of the paper is as follows. In sections 2 and 3 we review 
the von Neumann and Liiders versions of the projection postulate. In section 4 we 
present a reasonably self-contained account of the basic principles of stochastic state 
reduction, along with a brief synopsis of the relevant mathematical tools of stochastic 
analysis. For the purposes of illustration we consider primarily the case of an energy- 
based reduction model, for which the associated dynamics are given by the stochastic 
differential equation ([H]), though most of the relevant mathematical and physical ideas 
can be readily generalised to the class of reduction processes noted above, based on a 
compatible family of observables that commute with the Hamiltonian. 

Throughout the discussion we emphasise the role of martingale methods as an aid 
to the advancement of our understanding of quantum phenomena. In particular we 
show that the expectation of the Hamiltonian, which fluctuates during the course of the 
reduction process, is a martingale, and its variance is a super martingale. The martingale 
property satisfied by the expectation of the Hamiltonian can be viewed as a kind of weak 
conservation law for the energy, generalising the Ehrenfest relation. In section 5 we use 
the Doob-Kolmogorov maximal inequalities to obtain a set of upper and lower bounds on 
the fluctuations of the energy during the reduction process, and show that the magnitude 
of a typical fluctuation is roughly of the order of the initial energy uncertainty. We also 
demonstrate that the process followed by the squared energy uncertainty is given by the 
conditional variance of the terminal value of the energy. 

In sections 6 and 7 we derive the projection postulate, and prove that the collapse 
to the general Liiders state occurs with the appropriate probability, whether or not the 
initial state is pure and whether or not the eigenvalue spectrum is degenerate. Then in 
section 8 we study the implied evolution of the density matrix for a given initial state, not 
necessarily pure, and derive an exact time-dependent solution for the Lindblad equation 
associated with the reduction. We conclude in section 9 by introducing a change-of- 
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measure technique to solve the stochastic differential equation for the dynamics of the 
state vector, thereby constructing an explicit unravelling of the Lindblad equation. 

2. Von Neumann Projection 

Let us consider a quantum system which, for mathematical simplicity, we shall assume 
is characterised by a finite dimensional Hilbert space 7i of dimension N. Suppose that F 
is an observable for which the corresponding Hermitian operator acting on Ti is denoted 
F. Our interest here is in providing a clearer understanding of what happens if the 
observable F is measured when the system is in a given pure state, corresponding to a 
ray through the origin in 7i. 

First we shall examine the more straightforward case where F has a nondegenerate 
spectrum, and the eigenvalues of F are given by the numbers f n for n — 1, 2, • • • , N, with 
the property that /„ ^ f m for n ^ m. The eigenstate corresponding to the eigenvalue 
/„ will be denoted |/ n ), so that F\f n ) = f n \f n ). 

Let us write \ipo) for a representative state vector for the given initial pure state. 
We shall assume that \ip ) is normalised to unity, so (ipo\ipo) — 1- We shall also assume 

tnat \fn\jm) — Onm- 

According to the projection postulate, when a measurement of F is made, the state 
vector undergoes a transition \ip ) \— > \f n ) to one of the eigenstates of F. This occurs 
for a specified value of n with the probability 

^n=|(/n|^0)| 2 , (1) 

and the result of the measurement in that case is the eigenvalue f n . The associated 
transition \ipo) i— > \f n ) is called the state reduction or 'collapse of the wave function' 
arising from the measurement of F. 

A more precise way of stating this is that when the observable F is measured, 
the initial pure state |f/>o) transforms to a 'mixture', i.e., a random state |f) with the 
property that |f) is given by the eigenstate \f n ) with the probability n n . The fact that 
|f) is a mixture reflects our ignorance of what the result of the measurement process 
will be. When the observable F is measured, we can be confident that the result is one 
of the values f n , and that the new state is the eigenstate represented by \f n ), but we 
cannot say in advance which one it will be. This is what is meant by saying that the 
result of the measurement is random. 

The density matrix p associated with the random state |f) and the probability 
distribution 7r n is given by the expectation of the random projection operator |f)(f| 
associated with |f). In other words, we have 

p = E[|f)(f|] 

N 

= 5>rJ/"X/n|- (2) 

n=0 

Here E denotes expectation with respect to the distribution ir n . The mixture |f ) typically 
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carries more information than the density matrix p alone, because if one is given p, then 



in general there are many different mixtures that correspond to it ||12|| . 

The significance of the density matrix p is that if G is any other observable, not 
necessarily compatible with F, and we measure G after we measure F, then the expected 
value of G is 

E[(f|G|f>]= J> B (/„|6|/ B ) 

n 

= Tl J2*n\fn)(fn\G (3) 

n 

= TipG, 

where Tr denotes the trace operation. 

An alternative way of facilitating the description of the density matrix associated 
with the measurement outcome for an observable F is to introduce the projection 
operator P n = \f n )(fn\ associated with the eigenvalue f n . When F is measured, the 
state |^o) is transformed to 7r n P n \ipo) if the measurement outcome is known to be 
f n , for which the corresponding probability is 7r n . An analogous transformation holds 
for the corresponding density matrix p — l^oX^ol which transforms according to the 
scheme 

PO >-* Poo = Kn lp nf>0 P n 

= \fn}(U (4) 

Here we use the notation po to signify the density matrix before the measurement, and 
Poo to signify the density matrix after the measurement. 

More generally, if F is measured but the outcome is not known, then po transforms 
according to the scheme 

N 
/ J PnP0 P n 



PO ^ P 



oo 



n=l 

N 

= ^|/n)(/n|^0)(^0|/n)(/n| (5) 

n=l 

N 

n=\ 

The density matrix itself is often referred to as representing the 'state' of a quantum 
system. This is because it contains all the information required to calculate ensemble 
expectations and probabilities for the measurement outcomes of quantum observables, 
conditional on the present state of knowledge of the observer. 

In the examples above, the ensemble interpretation is as follows. We prepare a large 
number of independent identical quantum systems each in the state IV'o); then measure 
F . Now there are two possibilities. In the first case, we only keep those systems for 
which the result of the measurement of F was the value f n . The new ensemble then 
consists of a large number of independent systems each of which is in the pure state 
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\f n ). The corresponding density matrix is \f n ){fn\- I n the second case, we keep all the 
systems after F has been measured. The resulting ensemble is therefore a mixture, and 
for each value of n, a given system is in the pure state \f n ) with probability ir n . The 
corresponding density matrix is then given by (|5|). 

3. Liiders' Postulate 

Let us now turn to the case of an observable with a degenerate spectrum. In this case, 
we shall write \f n j) for an orthogonal basis of distinct eigenstates of F sharing the same 
eigenvalue f n . Here n = 1, 2, • • • , D, where D is the number of distinct energy levels, 
and j = 1, 2, • • • , d n , where d n is the dimensionality of the subspace 7i n of 7i spanned 
by the eigenstates with eigenvalue f n . For convenience we normalise \f n ,j) such that 

(fn,j\fm,k) = S nm 5jk. (6) 

Then the projection operator P n onto the subspace spanned by states for which F = f n 
is given by 

Pn = J2\f^)(fn,\- (7) 

We note that P n is independent of the specific choice of basis made for the designated 
subspace, and that P n P m = P n 5 nm and P n \f m ,k) = $nm\fm,k)- 

With these preliminaries in mind, suppose the observable F is measured when 
the system is in the pure state |^o) 5 and the result of the measurement is one of 
the degenerate eigenvalues f n . In this case, it is perhaps less obvious a priori what 
the correct probability is for the outcome, or indeed what becomes of the state once 
the measurement result is known. A more refined version of the projection postulate 
is required to deal with this situation, due to Liiders 0, according to which the 
measurement outcome probability is 

Prob[F = f n ] = (ij \P n \tp ), (8) 

and the associated state reduction is given by 

|V0> ^ K- 1/2 Pn\*Po), (9) 

where 

7T n = (lff \P n \lpo) 

d n 

3=1 

Thus, of all the possible eigenstates with eigenvalue f n , a single choice is made, given 
by the projection from the initial state vector onto the relevant subspace. As in the 
nondegenerate case, the measurement outcome can be described by a mixture, where 

-I Icy ^ 

the random state |f) is given by the normalised Liiders state \f n ) = n n P n \ipo) with 
probability 7T n (n = 1, 2, • • • , D). 
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The validity of the Liiders postulate can, in principle, be tested by a succession 
of measurements of the energy of a system, followed by the measurement of another 
observable incompatible with the energy. Consider for example the system consisting 
of a pair of noninteracting spin-| particles in an external magnetic field aligned along 
the z-axis, for which the energy eigenstates are the spin-0 singlet and spin-1 triplet. 
The corresponding eigenvalues are degenerate and are given by E\ — — 1, E^ — -E3 = 0, 
and £4 = 1, in suitable units. We choose the initial state such that the Liiders state 
associated with the degenerate energy eigenvalue is given by the spin-0 singlet. Then, for 
an ensemble of identically prepared systems, we measure the energy, and discard those 
systems for which the outcome is given by one of the eigenvalues ±1. The remaining 
systems, according to Liiders, are in the spin-0 singlet state, whereas according to von 
Neumann fl[], these states would in general be in superpositions of the singlet and S z = 
triplet states, the precise details of which depend on the nature of the measurement 
apparatus. In the present set-up, one can measure the total spin operator to determine 
the outcome of the initial energy measurement. This is because the result of the total 
spin gives the eigenvalue if and only if the system is in the singlet state. 

An important feature of the Liiders postulate is the inherent 'instability' of the 
reduction process implied for certain types of measurements. That is, in the case 
of an observable with a degenerate eigenvalue f n , the projection is onto a single 
state ii n P n \ipo)] whereas if the observable is perturbed even slightly, breaking the 
degeneracy and producing, say, two distinct but close eigenvalues f n and f n ', then the 
reduction process bifurcates, leading to one or the other of two orthogonal, and hence 
maximally separated, eigenstates |/ n ) and \f n >). 

Thus we are led to consider that if the phenomenon of state reduction itself arises 
as a consequence of a dynamical process, then this process must have sufficiently 
special properties to ensure that under a smooth deformation of the parameters 
characterising the observable being measured, the resulting dynamics exhibits the 
required discontinuous behaviour and produces the corresponding bifurcation in the 
postulated transition probabilities. In what follows, we shall demonstrate that the 
standard stochastic models for quantum state reduction exhibit the required special 
properties, including the relevant bifurcation phenomena. In short, we can derive the 
projection postulate from a dynamical model, including the specifics of the Liiders 
postulate and its consequences in the case of a degenerate spectrum and a mixed initial 
state. 

4. State Vector Reduction 

In this section we consider in some detail the case of the standard energy-based stochastic 
extension of the Schrodinger equation, for which, if we set H = 1, the dynamics are given 
by the following stochastic differential equation on 7i: 

d|Vt) = - iH\if> t )dt - \o-\H - H t ) 2 \^ t )dt + \a{H - H t M)dW t . (11) 
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The properties of dynamical processes of this type have been investigated by a number 
of authors ||, ||, [7|, §, [| . In particular, as we shall demonstrate with various examples, 
many of the familiar probabilistic features of standard quantum mechanics, including 
the Born rules, can be deduced from (PI), or suitable generalisations thereof. Here, \ip t ) 
denotes the random state vector at time t, for which an initial state |^o) is prescribed. 
For the moment we consider the case when |^ ) is known, and later we turn to the case 
where |^o) is random. For the expectation of the Hamiltonian operator H in the state 
\ipt) we write 

(Aim) 

mm) 

from which it follows that H t itself can be interpreted as a random process, which we 
shall call the energy process. Strictly speaking, an expression such as (H — H t ) 2 \ip t ) in 
(|TT|) should be written (H — H t l) 2 \ip t ), where 1 is the identity operator, but there will 
be no ambiguity if we use the more compact notation. 

The energy-based stochastic extension of the Schrodinger equation (|TT| ) is of great 



interest because it represents perhaps the simplest plausible model for the collapse of 
the wave function, and as such exhibits many remarkable features, both physically and 
mathematically. It will be useful if we begin our analysis with a brief overview of the 
probabilistic framework implicit in the characterisation of flTTl). 

The stochastic differential equation for the process \tp t ) is to be understood as 
defined on a fixed probability space (Q, J 7 , P) equipped with a filtration JF 4 , with respect 
to which W t is a standard Wiener process (Brownian motion). Here Q is the sample 
space, T is a a- field on Q, and P is a probability measure on T . 

The filtration Tt represents the information available at time t, where t G [0, oo). 
More precisely, a filtration of T is a collection !F t (0 < t < oo) of a-subfields of J 7 
with the property that s < t implies T s C Tt- A function X : Q i— > K. is said to be 
measurable with respect to T if for each x G K. the set consisting of all u G Q satisfying 
X(u) < x is an element of T . This assures that Prob[X < x] exists with respect to 
the given measure P on J 7 , and we say that X is a random variable on the probability 
space (fi, J 7 , P). Then by a random process we mean a parameterised family of random 
variables X t (0 < t < oo) on (Q, J 7 , P). With a slight abuse of notation we let X t stand 
both for the entire process X t (0 < t < oo), as well as the random variable X t for 
some given value of t; usually it will be evident from context which meaning is intended. 
Likewise T t may denote the entire filtration Tt (0 < t < oo), or the a-subfield of J 7 
corresponding to the information set at time t. If a random process X t is such that for 
each value oft the corresponding random variable X t is jF r measurable, then we say that 
the process X t is adapted to the filtration T t - The idea of 'adaptedness' is important 
because it is though this device that a notion of causality is introduced for the class of 
process we consider. 

It should be emphasised that the probabilistic concepts outlined here and in what 
follows are introduced not merely for the sake of mathematical clarity (although this 
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is in itself a desirable feature), but also because it makes possible the use of various 
powerful analytical tools, examples of which we shall discuss shortly. 

The concept of conditional expectation plays a particularly important role in the 
theory of quantum state reduction, and hence it will be helpful if we elaborate slightly 
on the idea here. For any random variable X on (fl, JF, P) we define its expectation 



E[X] = / X(uj)dP(uj) (13) 

by use of the Lebesgue integral. One can think of E[X] as the ensemble average of 
X. Then if £ is a u-subfield of J-, the random variable Y is said to be (a version 
of) the conditional expectation of X with respect to £ if Y is ^-measurable and if 
E[X1 A ] = E[Y1 A ] for all sets A e £. In that case we write Y = E[X\£]. Here l A 
denotes the indicator function for the set A, so 1,4(0;) = 1 for u E A and 1,4(0;) = for 
uj (£ A. By convention we write E f [X] = E[X|jF t ] for the conditional expectation of X 
given information up to time t. One can think of E*[X] as the ensemble average of X 
conditional on the history of events up to time t being specified. 

A useful result that follows on from these definitions is the so-called 'tower 
property' of conditional expectation, which says that if "D is a a-subfield of £ then 
E[E[X|£]|X>] = E[X|T?]. If we set V = (0,0), the smallest a-subfield of T, then 
E[X|X>] = E[X]. This is because the only D-measurable random variables in that case 
are the constants: if X(u>) is constant on Q, then for any given x we have X(u>) < x 
either for all uj G Vt or for no w G fl; conversely, if X(uj) is not constant, then we can 
find a value of x and two points w^^Gfl such that X(ui) < x and X(o^) > x, which 
shows that X(o;) is not P-measurable — that is to say, the set {o; : X(o;) < x} is not 
an element of T>. It follows then from the tower property that E [E[X|£]] = E[X], the 
so-called law of total probability. In the case of a filtration Tt (0 < t < 00), if we 
take £ = J-'t and V = JF S , then the tower property reads E s [E t [X]] = E S [X] for s < t; 
whereas the law of total probability implies that E [E 4 [X]] = E[X] for all t > 0. 

Now suppose X t (0 < t < 00) is an adapted process on a probability space 
(il,^ 7 , P) with filtration Tt (0 < t < 00). Then we say X t is a martingale if the 
following two conditions hold: E[|X t |] < 00 for all < t < 00, and E s [X t ] = X s for all 
< s < t < 00. If instead of the latter condition X t satisfies E s [Xj] > X s then we say 
X t is a submartingale, and if E S [X 4 ] < X s we say X t is a supermartingale. 

Returning now to our investigation of the process (|TT|), we note that the parameter 
a governs the rate at which the state vector reduction proceeds for a given level of initial 
uncertainty in the energy. The units of a are 

a ~ [energy] -1 [time] -1 ' 2 . (14) 

The characteristic time-scale tr associated with the collapse of the wave function is 






a 2 V n 



where Vq is the square of the initial energy uncertainty AH. 
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In what follows we shall make no specific assumptions about the value of a . 

-| Icy 

Nevertheless, we note, as is discussed in J7|, |[, that if a ~ M p in microscopic units 
with h = c = 1, where M p is the Planck mass, then the 'large numbers' cancel out and 
we are left with a typical reduction time-scale of 

m ~ {-^r) b. (16) 

This expression is interesting as a candidate for t r inasmuch as it relates energy spreads 
typical of atomic and nuclear phenomena to time-scales that are accessible in the 
laboratory. 

The factor of ~ appearing in front of a in (JTTJ) is for convenience, and ensures 
consistency with the notation of J7|, §, ||. 

It follows from the Ito rules (dt) 2 = 0, dtdW t = 0, and (dW t ) 2 = dt, as well as the 
special form of the nonlinear terms appearing in (|Tl~D , that the norm of the state l^) is 
preserved under the evolution (|TT|) . This can be seen as follows. The Ito product rule 
states that if X t and Y t are Ito processes then d(X t Y t ) = Y t dX t + X t dY t + dX t dY t . As 
a consequence, we have 

d«tk|V*» = (d(^l)l^) + mWt)) + (d(^|)(d|^». (17) 

Now the Hermitian conjugate of (|TT| ) is 

d(Vt| = i(A\Hdt - y(^t\(H - H t fdt + \a(il> t \(H ~ H t )dW t (18) 

Therefore, by use of the Ito rules, we obtain 

(d(^|)|^> = (iH t - yv t ) (Vtl^)dt, (19) 

and its conjugate, 

(^|(d|V t )) = {-iH t - \a 2 V t ) (VtlVt)dt, (20) 

together with 

(d(^|)(d|^» = iaV t (^|^)dt, (21) 

where V t is given by formula (^3) below. It follows then from ( |T7D that d(( , i |'?/' 4 )) = 0, 
as desired. This result is useful in calculations because we can assume the initial norm 
to be unity, without loss of generality, and thus {i>t\i>t) = 1 for all t. 

An analogous calculation shows that the energy process H t defined in (O) satisfies 

dH t = <rV t dW t , (22) 

where V t is the process for the variance (squared uncertainty) of H in the state \ipt), 
given by 

Vt ~ mh) ■ (23) 

The variance process for the Hamiltonian has the property that Vt = at time t if and 
only if \ip t ) is an energy eigenstate at that time. As a consequence of (|2"2]), we can write 

H t = H + a [ V u dW u , (24) 



o 
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where H Q is the initial expectation value for the energy. Now it is a general property of 
the stochastic integral that for any ^-adapted integrand A t satisfying E J A\du 
we have 



< oo 



E a 



A u dW u 
o 



A u dW u (s<t). (25) 

'o 

The variance process V t is bounded by \{E + — E_) 2 where E + and E- are the largest 

and the smallest energy eigenvalues, respectively, which implies that f Q V 2 du < oo. 

Furthermore, we note that \H t \ is bounded by max(|.E + |, |-EL|). It follows that H t is a 

martingale: 

E s [H t ] = H s , (s<t). (26) 

The martingale condition is the stochastic analogue of a conservation law, and thus ( p6|) 
can be interpreted as a weak conservation law for the energy. We recall that, for ordinary 
quantum-mechanical evolution in the case of a time-independent Hamiltonian, the 
Schrodinger equation dt\ipt) — —iH\ip t ) ensures that the expectation of the Hamiltonian 
Ht = (ipt\H\ipt) / (iptlipt) is conserved along the Schrodinger trajectories. In the case of 
the stochastic extension of the Schrodinger equation we have instead the martingale 
relation (^) which ensures that the ensemble average of the energy is conserved. 

Because ( fffi) plays a pivotal role in understanding the nature of the reduction 
process, we shall sometimes refer to the system of stochastic dynamics described by 



( p|) as a martingale model for quantum state reduction. 

We note, more generally, that if the operator G = g{H) is given by a function of 
H, then the process 

(Wt> 
is also weakly conserved, i.e., E g [G t ] = G s for s < t. Thus, for example, if g{x) = x n 
and we introduce the notation 

Hi* - <™ (28, 

mm) 

for the n-th moment of the energy, then 

dH[ n) = a (H\ n+l) - H t Ht ] ) dW t , (29) 

where H t — H\ . 

With these formulae in mind, let us consider now the dynamics of the variance 
process V t . Writing V t = HJ; ' — (H t ) 2 it follows, according to Ito's lemma, that 

dV t = dHi 2) - 2H t dH t - (dH t f. (30) 

By use of the Ito rules together with ( |2"9"D we then deduce that 

dV t = -a 2 V t 2 dt + o-fadWt, (31) 



where 



a - mi^iptA (32) 



Martingale Models for Quantum State Reduction 



11 



is the skewness of the energy distribution at time t, i.e., the third central moment of 



the Hamiltonian. More specifically, we have (3 t — H\ — ?>H t H). + 2 (if, 



(3) 



r(2) 



Integrating equation fl3~l|) for the dynamics of the variance we obtain 

V = V - a 2 [ V 2 du + a [ (3 u dW u , (33) 

Jo Jo 

from which it follows at once, by use of (^5|), that 



and thus 



E s [V t } = V s -a 2 E s 



E s [V t ] < V s 



VMu 



(34) 



(35) 



for s < t, which shows that V t is a supermartingale, i.e., a process that on average 
decreases. 

In particular, if we write Vt = E[Vj] for the ensemble average of Vt, then it follows 
as a special case of ( p4|) that 

,-t 



V = Vn 



a 2 E 



v 2 du 



Uo 



(36) 



(37) 



Differentiating this expression with respect to t we obtain 

where the process i] t defined by i] t = E[(V t — V t ) 2 ]/V 2 is nonnegative. Therefore, by 
integration of (|37|) , we obtain 

V 



V 



l + a 2 K,(t + 6)' 
where ^t = J Vsds. As a consequence we deduce that 

V 



V,< 



which shows that 



l + a 2 V r 



lim V = 0. 



t— KXD 



(38) 



(39) 



(40) 



Alternatively, if we introduce the 'localisation' process A t = Vf 1 then it follows from 
(|37l) that d t A t > a 2 , which shows that A t increases without bound [14], |1^]. In (p^) 
we see an example of the role of 77? = (cr 2 Vo) _1 as the characteristic time-scale of the 
reduction process. Since V t is nonnegative, it follows that 

lim V t = (41) 

t— >oo 

almost surely. The dynamical process (|ll|) therefore induces a collapse of the wave 
function, for any choice of the initial state \ipo), to an eigenstate of the Hamiltonian. 
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The martingale property ( |2q ) satisfied by the energy process H t implies in the limit 
t — > oo, that Efi^oo] = H . However, the terminal value H^ of the energy process, the 
existence of which we shall establish shortly, is necessarily one of the energy eigenvalues, 
from which it follows that 



Hn 



2 , ^nEr, 



(42) 



where 7r n is the probability of reaching the eigenstate \n) starting from the given initial 
state. Therefore, the ensemble average of the measured value of the energy equals the 
expectation value of the energy in the initial state, as it should. 

The importance of this conclusion is that whereas in quantum measurement theory 
it is essentially an assumption that the 'expectation value' of an observable in a given 
state is the ensemble average for the result of a measurement of the observable, in a 
martingale model one can prove that the asymptotic ensemble average agrees with the 
expectation value, hence justifying the conventional interpretation of this quantity. In 
particular, using (|26| ) we can write 

H t = E t [H OQ ] (43) 

which shows that the quantum expectation value H t of the observable H at time t is 
always the best 'forecast', based on information currently available, for the outcome of 
a measurement of H. 

A similar result holds for the dispersion of the measured values of the energy. This 
can be established by use of the Ito isometry. If A t and B t are JF r adapted real processes 
that are square-integrable in the sense that E[f Q A 2 ds] < oo and E[f Q B 2 ds] < oo, then 
the Ito isometry states that 



E 



AAW K 



BAW* 



E 



A.BAs 



uo 



(44) 



It follows therefore from (|24j ) that 

E [{H t - H ) 2 } 



a 2 E 



a 2 E 



(/'w,) 2 


Jo 





(45) 



by virtue of the Ito isometry. By use of expression (^3|) for V t we then deduce that 

E[(H t -H ) 2 } =V -E[V t }. (46) 

Taking the limit t — > oo and using the fact that lim^oo V t = 0, we get 

Varttfoo] =E[(H 00 -E[H 00 }) 2 ] 

= V , (47) 

which demonstrates that the variance of the measured energy is in agreement with the 
squared energy uncertainty in the initial state. 
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During the course of the reduction process, the energy H t of the system can, in 
principle, deviate far from its initial value H Q , subject to the condition that it stays in 
the range H t G [E-,E+], where EL and E + are the lowest and highest energy levels. 
Nevertheless, we can show that on average H t will not deviate too much from H : an 
upper bound can be set on the maximum fluctuation experienced by the energy, on 
average, as the reduction proceeds. This bound is given by 2AH, twice the initial 
energy uncertainty. 

The proof of this result makes use of the Doob-Kolmogorov maximal inequalities 
(see, e.g., fL6|, theorem 6.10, or flTfl, theorem 1.7, p. 54). These inequalities state that 
if M t is a right-continuous martingale or positive submartingale, and E[|Mt| p ] < oo for 
some p > 1, then 

p 



E 



sup 

0<t<T 



\M t \ p 



< 



V 



p — 1 



E[|M r | 



and 



Prob 



sup \M t \ > K 
0<t<T 



< — ENMtH 



(p>l), 



(p>l) 



(4* 



(49) 



for any constant k > 0. 

In the present context, we are especially interested in the inequality obtained in 
the case p = 2, for which we have the relation 



E 



sup M 2 t 



0<t<T 



< 



4E [M%] 



(50) 



which is known as Doob's L 2 -inequality. Now, setting M t = H t — Ho and using equation 
(filf) we obtain 



E 



sup (H t - H ) 

0<t<T 



<4(V -V T ). 



In particular, taking the limit T — > oo, it follows from (£|l|) that 



E 



sup (H t 

0<t<oo 



Hn 



<4^ , 



(51) 



(52) 



which shows that, on average, the energy stays within two standard deviations of its 
original value. 

This result is consistent with the intuition often arising in physical arguments to the 
effect that when a system is in a state of uncertain energy, then the energy fluctuates, 
with a typical fluctuation roughly of the magnitude ~ AH. There is no quantum- 
mechanical principle which states that such fluctuations actually occur, but one can 
see that in a martingale model there may indeed be a natural basis for inferring the 
existence of fluctuations of the required magnitude. We note that the bound implied by 
the inequality ( |52|) is independent of the choice of a, which shows that it is valid also 
for relatively stable, long-lived states, i.e., those for which cr 2 Vo is small. 

^From ( pE9"D we can determine an upper bound on the probability that the magnitude 
of the energy fluctuation will exceed any designated threshold during the reduction 
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process. Specifically, if we set p 
obtain 



2 and k = X^/Vq, then taking the limit T — > oo we 



Prob 



sup (H t 

0<«oo 



Ho) 2 > AVq 



1 



(53) 



A related bound for the variance process Vt can be obtained by use of Doob's 
maximal inequality for positive supermartingales (see, e.g., [17 1, p. 58). This relation 
states that, if X t is a right-continuous positive super martingale, then for any constant 
k > 0, we have 



Prob 



sup X t > k 

0<«oo 



< ^E[X ] 



(54) 



In the case of the variance process Vt, which as we have shown is a positive 
supermartingale, if we set k = \ 2 V , then fl54]) becomes 

1 



Prob 



sup V t > A V 

0<«oo 



< 



A 



2' 



(55) 



This relation shows that, during the reduction process, although the energy variance 
can increase owing to random fluctuations, there is a bound on the probability that the 
energy uncertainty ever reaches A times the ensemble average of the initial uncertainty 
for any given value of A, and this bound is given by A~ 2 . 

Let us return now to the asymptotic relation Efii/oo] = H and ask whether the 
terminal value H^ of the energy process actually exists as a random variable. To 
prove that it does, we make use of the martingale convergence theorem, which in a form 
sufficient for our purpose states that if a continuous martingale M t satisfies E [|M 4 | P ] < k 
for some p > 1 and k < oo, and for all t G [0, oo), then there exists a random variable M^ 
satisfying E [|M 00 | P ] < k and M t = 'E t [M 0O \, with the properties that lim^oo M t = M^ 
almost surely and that lim^oo E [\M t — Moo| p ] = 0. 

In the present context, by setting M t = H t — H , we thus deduce the existence of 
an asymptotic random variable H^ with the property that H t converges to H^ almost 
surely, and H t = Et [#«,]. 

We conclude this section by generalising (|47|) to demonstrate that the energy 
variance process V t defined by (|23[) has the natural interpretation 

V = Var 4 [H^] . (56) 

That is to say, V t is given by the conditional variance of the terminal value of the energy, 
given information up to time t. To establish this relation we proceed as follows. 

For any random variable X on the probability space (Q, T , P) we define the 
conditional variance Var[X|£] with respect to the a-subfield £ C T by 

Var[X|£]=E[(X-E[X|£]) 2 |£]. (57) 

It follows as an application of the law of total probability that 

Var[X] = E [Var[X|£]] + Var [E[X|£]] , (58) 
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the so-called conditional variance formula. Thus for example if J-" t (0 < t < oo) is a 
filtration of (fl, J 7 , P) and we write Var t [X] = Var[Jf|jF t ], then 

Var t [X]=E t [(X-E i [X]) 2 ] (59) 

and for the conditional variance formula we have 

Var[X] = E [Var t [X]} + Var [E t [X]} . (60) 

In the problem at hand, we note that in the limit t — > oo formula ( |33D for the 
variance process V t takes the form 

POD 

(61) 



OO /"OO 

2 / t^2. 



Vb + o- / /3 u dW n = a z / V u z du. 
Jo Jo 

Therefore, taking the conditional expectation of each side of this relation and using 

formula (|5|) we deduce that 



V Q + cr P u dW u = a 2 E t 
Jo 



VMu 



(62) 



Substituting this relation into (1331) then gives us 



Vt = v 2 E t 



a 2 E t 



o- 2 E t 



V 2 du 



VMu 



V u dW u 



rr i V u z du 

'0 



(63) 



= E t [{H^ - H t ) 2 ] 
= Var t [Hgo] , 

as desired. We note that in the next to last step here we have used (|24l) together with 
the conditional Ito isometry 



E t 



A u dW u 



E t 



A 2 u du 



(64) 



< oo. 



valid for any adapted integrand A u satisfying E f Q A^du 

A positive supermartingale with the property that its expectation goes to zero 
asymptotically is called a potential [Q. The analysis above shows that the variance 
process associated with quantum state reduction satisfies these conditions and admits 
a Doob-Meyer decomposition of the form 

V t = E t [Z^] - Z t , (65) 

where 



a 



V 2 du 



(66) 



o 



is an increasing process. 
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6. Reduction Probability 

The probability ix n of reduction to a specific energy level E n , under the dynamics 
governed by the stochastic differential equation (|TTD , can be determined as follows. 
The method we use is essentially that of || ||. 

First we observe that for any operator G acting on 7i, the process Gt for the 
expectation value of G in the state \ip t ) satisfies 

dG t = -i(^ t \[G,H}\^ t )dt+y^ t \ (HGH-\{H\G}) [Vt>dt 

+ la(ipt\{(G - G t ), (H - H t )}\ij t )dW t . (67) 

Here [X, Y] = XY - YX and {X, Y} = XY + YX denote the commutator and the 
anticommutator, respectively. 

The drift term in ( |6^) consists of two parts: the first is the familiar Ehrenfest term 
involving the commutator with the Hamiltonian; the second is a term of the Lindblad 
type [H, [G, H]] arising as a consequence of the diffusive dynamics of the state vector. 
The volatility term in fl6"T|), i.e., the coefficient of dWt, is given by the covariance of G 
and H in the state \ipt)- If G and H commute, then the drift vanishes, and fl6"T|) reduces 
to 

dG t = a (jtftlGHtyt) - G t H^j dW t , (68) 

from which it follows that the process G t is a martingale || . This is consistent with our 
earlier observation that H t is itself a martingale, and that the process G t corresponding 
to any function of the form G = g{H) is also a martingale. 

Now let us consider the projection operator P n for the subspace Ti n of Tt spanned 
by the energy eigenstates with energy E n . In the case of a nondegenerate eigenvalue, 
we have P n = \n)(n\. On the other hand, if E n is a degenerate eigenvalue, then 

Pn = Y,\ n ^){^3l (69) 

i=i 

as in (0), where d n is the dimension of the subspace H. n and \n,j) (j = 1, 2, • • • , d n ) is 
an orthonormal basis for 7i n . Clearly P n commutes with the Hamiltonian H for any 
value of n. Furthermore, the relations HP n = P n H = E n P n and H = J2 n EnPn are 
equivalent on account of the resolution of identity 

D 

E^ = 1 > ( 7 °) 

n=l 

where D is the number of distinct energy eigenvalues. 
Now let us write 

mm) 

for the expectation of the projection operator P n in the state \ip t )- Because P n commutes 
with the Hamiltonian, we deduce that the process P nt is a martingale for each value of 
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n. We note that J2 n P nt = 1 and J2 n E n Pnt = Ht- In particular, by setting G t = P n t, 
one infers from (|68| ) that 

dP nt = aP nt (E n - H t )dW t . (72) 

This stochastic differential equation implies that P nt will continue to fluctuate as long 
as H t ^ E n and P nt ^ 0. The solution of ( [72] ) is given by P nt = P n oM nt , where 

M nt = exp fa [ {E n - H s )dW s - \a 2 f {E n - H s ) 2 ds\ , (73) 

and P n o is the initial expectation value of the projection operator P n , This follows 
from the fact that for any bounded ^-j-adapted process at the solution of the stochastic 
differential equation dX t = a t X t dW t (X > 0) is 

X t = X exp (J a s dW s - \ j a 2 s ds\ , (74) 

which one can verify by an application of Ito's lemma. Because P nt is a martingale, it 
follows that 

E [P nQO ] = P n0 . (75) 

Here E[P noo ] is the ensemble average of the expectation value of the projection operator 
P n at the terminal energy eigenstate of the reduction process. Because P noo takes the 
value one if the terminal energy has eigenvalue E n and takes the value zero otherwise, 
it follows that E[P noo ] is the probability of reaching a state with energy E n , i.e., 

E [P noo ] = Tl n . (76) 

With these observations at hand, we are now in a position to interpret the 
asymptotic martingale relation (|75|) . If E n is a nondegenerate eigenvalue, then P n0 



is the usual expression for the Dirac transition probability from the initial state \ipo) to 
the eigenstate \n), given by 

(ip \n){n\ip ) 

{tpo\ipo){n\n) 
Thus we conclude, in the case of a nondegenerate Hamiltonian, that the martingale 
model for quantum state reduction allows one to deduce the correct transition 
probabilities. 

In the case of a degenerate eigenstate, the probability 7i n can also be interpreted in 
terms of a Dirac transition probability. In particular, whether or not the spectrum of 
the Hamiltonian is degenerate, we can write 



■n 



«Vo|Pn|V>0» 2 _ 
/>0> 

(7f 



(lpo\lpo) (lpo\Pn\lpo) 
(P n j'o\^o)(^o\Pn^o) 
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_ (P n 1po\lpo)(lpo\Pn1po) 

(MMMPZliM 

_ (Pni>o\i>o)(4>o\Pn'>J>o) 

(lpo\lp )(Pni'o\Pn'4'o)' 

where the Liiders state \P n ip) is defined by 

\P n ^)=Pn\^)- (79) 

Therefore, by virtue of (]75|) and (|78|) , we see that the probability of obtaining the 
eigenvalue E n is given by the Dirac transition probability from the given initial state 
\ipo) to the Liiders state \P n iJjo). 

The interesting point here is that, once again, while this is an assumption in 
standard quantum theory, it arises as a deduction in the martingale model for quantum 
state reduction. 

In fact, we can demonstrate, in the case of a degenerate eigenvalue, that the 
reduction necessarily results in the Liiders state if the corresponding eigenvalue is 
obtained. This can be seen as follows. 

For each value of n — 1, 2, ■ ■ • , D, such that P n \ipo) ¥" let us write 

|n, 1) = 5^ (80) 

(^olWo) 1 / 2 

for a basis vector corresponding to the normalised Liiders state for that projection 
operator, and let \n,j), j ^ 1, be an associated basis for the states orthogonal to \n, 1) 
that lie in the subspace Tt n spanned by eigenstates of energy E n . The operator 

ftn = ^2\n,j){n,j\ 

= P n -\n,l)(n,l\ (81) 

thus projects onto the subspace of H n consisting of vectors orthogonal to the Liiders 
state for that value of n. Evidently, we have 

rWo> = (82) 

which follows from fl8"0"D and the fact that II n -P n = Il ra . Since the projection operator Il n 
commutes with the Hamiltonian, the process 

(if} t \f[ n \ip t ) 

mm) 

is a martingale, the initial value of which is n n0 = on account of the relation 
flB2|) . Therefore, by virtue of the martingale relation E [n noo ] = n n0 , we deduce that 
E[II noo ] = 0. Now, Il noo is a nonnegative random variable. Therefore, if E[n noo ] = 
then Il noo = almost surely. It follows that the terminal state must be orthogonal to 
the subspace of 7i n spanned by states with energy E n that are orthogonal to the Liiders 
state. As a consequence, we see that if reduction occurs to a state of energy E n , then 
that state must be the Liiders state corresponding to that eigenvalue. 
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In fact, we deduce a stronger result: namely, that the stochastic motion of the 
state vector, during the course of the reduction process, is necessarily confined to the 
D-dimensional subspace ofH spanned by the Liiders states P n \ipo) , for n = 1, 2, • • • , D, 
where D is the number of distinct energy eigenvalues and P n is the projection operator 
onto the subspace TC n ofTC spanned by eigenstates with eigenvalue E n . 

The proof of this theorem follows from the fact that, for each n, the process U nt is 
a martingale, and because il n0 = we have E[Il„ t ] = for all t > and thus U nt = 
for all t > 0. Therefore the state vector \i(j t ) always lies in the space spanned by the 
vectors P n \ipo) for n — 1, 2, • • • , D. 

A similar analysis is valid in the more general situation for which the dynamics of 
\i/) t ) are given by a stochastic differential equation of the form 

r 

d\i/H) = ~ iH\,k)dt - § J>a(£* - F <*t?\i>t)dt 

a=l 
r 

+ \Y,° a (F a -F a t)dW?. (84) 

a=l 

Here F a (a = 1, 2, • • ■ , r) represents a commuting family of observables, each of which 
also commutes with the Hamiltonian H, the o~ a are associated coupling constants, and we 
write F at for the expectation of F a in the state \ipt)- In this case W? denotes a standard 
r-dimensional Brownian motion, and the reduction proceeds to a common eigenstate of 
operators F a (a — 1, 2, • • • , r). Most of the results of this paper are applicable mutatis 
mutandis to this more general class of reduction process, though in what follows we shall, 
for simplicity, continue to confine the detailed discussion to the case of the energy-based 
reduction (|Tl"l). 

7. The Case of an Initially Mixed State 

Thus far we have considered the role of the Liiders postulate as it applies to an initially 
pure state \ipo), and we have demonstrated that the postulate follows directly as a 
consequence of the martingale model for quantum state reduction. The Liiders postulate 
is, however, applicable in a more general context as well: namely, when the initial state 
is specified as a mixture with density matrix po- In that case, when an observable F is 
measured, the associated state reduction is given by the Liiders rule 

Po " typT' (85) 

for the density matrix, if the measurement result is the eigenvalue f n , and this occurs 
with probability 

n n = Tr P n p . (86) 

Here, as before, P n denotes the projection operator onto the subspace H n of TC spanned 
by eigenstates with the eigenvalue f n , which may or may not be degenerate. 
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The interpretation of an expression involving density matrices, such as (pop, is best 
understood in terms of ensemble averages. Thus ( |S"5{ ) means that if initially p can 
be used to compute the expectation of any observable G, not necessarily compatible 
with F, then after F is measured, and if the result f n is observed, the density 
matrix P n p P n /TrP n p Q can be used to compute the expectation of G in a subsequent 
measurement. 

Additionally, given the initial density matrix p , if F is measured but no note is 
taken of the result, then the ensemble average for a subsequent measurement of the 
observable G is TrpooG, where 

Poo = Yl PnPoPn- (87) 

n 

It should be borne in mind that these expressions, while generally regarded as part 
of the standard apparatus of quantum theory, are not derivable from any of the more 
basic assumptions of quantum mechanics, and have to be regarded as constituting an 



additional postulate. See, e.g., |19] for an illuminating brief account of the status of the 
projection postulate in quantum mechanics, and its relation to state reduction. It is 
interesting to note that von Neumann, in his original splendid work on the subject |I|, 
apparently failed to offer a satisfactory expression for the density matrix in the case of 
the measurement of an observable with a degenerate spectrum, a deficiency only later 
rectified by Liiders and others || (cf. [p0| , section 9, and the remark attributed to A. S. 
Wightman on p. 550 of ref. |PH). 



The general Liiders rule (85) has the important property that, in the measurement 
of an observable with a degenerate spectrum, if the initial state is not pure, then the 
final state need not be pure, if the result of the measurement is one of the degenerate 
eigenvalues. 

Now let us see if we can gain a clearer understanding of the general Liiders 
formulae (|85|) , (|86D and ([37]) by consideration of the martingale model for quantum 
state reduction. In the theory of stochastic differential equations, it is acceptable that 
the initial value of the random process should itself be a random variable; thus it is 
merely a special case when |^o) in ( |TT| ) is known. The deterministic case corresponds 
to the situation where the initial density matrix is pure, i.e., of rank one. In the general 
case, where l^o) is random, i.e., given by a mixture, the corresponding initial density 
matrix p is the ensemble average 

p = E[|*oX*o|], (88) 

where \*&o) is a random initial state vector, which we assume to be normalised. Then 
for the final density matrix we have 

P n |tf >(*o|£n 



Poo = J^E 



7T„f*r 



(* |P n |* ) 

as a consequence of the reduction 



(89) 



i*o> - /T ,;' T o) u/2 , 00) 
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where 7T n (^/ ) is the conditional probability that the eigenvalue E n is obtained, given 
the random initial state |^o)- However, this conditional probability is given by 

7T„(* ) = (*o|4|*0>, (91) 

i.e., the Dirac transition probability to the random Liiders state determined by the 
random initial state \^o), i n accordance with (]75|). As a consequence we see that (^) 
simplifies to give 

Poo = J2 E [ P n\*0)(*0\Pn 
n 

= ^P n E[|* )(* |]P n (92) 

n 

— / J PnPoPn, 

n 

and thus we obtain (^3). Additionally we have 

Poo = / T^nPnooi (.""J 

n 

where 

PnfaPn m 

Pn °° TrP n { ' 

is the reduced or 'conditional' density matrix, given that the observer has knowledge of 
the result H = E n , and 

7r n =E[7T n (*o)] 

= E[(* |P n |* ) 

= E[TrP„|* )(*o|J (95) 

= TrP n E[|* )(* |] 
= TrP n p 

is the probability of this result. 

8. Dynamics of the Density Matrix 

To gain further insight into the case where the initial density matrix is not pure, we can 
make a computation of the dynamics for p t . This can be achieved by examination of 
the Lindblad equation associated with the stochastic differential equation (|TTD, which 
in this case turns out to be solvable. 

If we start with equation (^7]) for the dynamics of the expectation value of an 
arbitrary operator G in the state \ipt), and take the ensemble average 

E [G t ] = TrGp t , (96) 
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where p t = E[|^ t )(0 t |], we find that 
dE[G t ] = TrGdpt 

= - iTrp t [G, H]dt + \<y 2 Tip t (hGH - \H 2 G - ±GH 2 ) dt (97) 

= - iTrG[H, p t ]dt + \a 2 TrG (Hp t H - \p t H 2 - \H 2 p^ dt, 

where in the second equality we make use of the cyclic property of the trace operator. 
This relation has to hold for any observable G, from which it follows that 

dtpt = -i[H, Pt ] + \a 2 (Hp t H - \H 2 p t - \p t H 2 ^ , (98) 

where d t = d/dt. This is the general equation of the Lindblad type [p^ , p3| associated 
with the stochastic differential equation ([11]), for which \o~H is the corresponding 
Lindblad operator. 

Now we consider the problem of solving the Lindblad equation ( |9"%| ) subject to an 
arbitrary specification of the initial density matrix po- For convenience we switch to a 
Heisenberg representation in which the density matrix is defined by the operator 

f t ± e^pte'™. (99) 

This has the effect of removing the purely unitary part of the evolution. For the 
dynamics of f t we have 

d t r t = e iflt (d tPt )e- iAt + i[H,r t } : (100) 

and therefore 

d t f t = \a 2 [kf t E - \H 2 r t - \f t H 2 ^ . (101) 

Let us write P n for the projection operator onto the subspace 7i n . Then, because 
P n H = HP n = E n P n , if we multiply each side of equation (|101|) by P n on both the right 
and the left we obtain 

dtiPnhPn) = 0, (102) 

from which it follows that P n f t P n is a constant of the motion. In particular, we have 
P n roP n = P n f<x>Pn, and thus P n p P n = PnP^Pn, and therefore 

J2 PnPoPn = J2 P nPooPn- (103) 

n n 

Because the terminal state is necessarily a mixture of energy eigenstates we have 

Y^ PnPooPn = Poo, (104) 

n 

from which by use of (|103| ) we immediately infer the general form of the Liiders reduction 
postulate flS^j). 

To proceed further we define the operator matrix valued process 1Z nm {t) by 

n nm {t) 4 P n r t P m . (105) 
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For each of the values of n and m, 7£ nm (£) is a time-dependent Hermitian operator. Here, 
n, m — 1, 2, • • • , D, where D is the number of distinct energy levels. From equation ( |101| ) 
for the dynamics of f t we deduce, by use of the relation P n H = E n P n , that 

d t n nm {t) = -W(E n - E m ) 2 n nm {t). (106) 



The general solution of the ordinary differential equation ( |106| ) is given by 

TLnmit) = Knm(0) exp (-§cr 2 (£ n - E m )H) 

= P n f P m e^{-\a 2 {E n -E m )H). (107) 

On the other hand, by use of the resolution of the identity (|70|) it follows from ( |105| ) 
that 

r t — / T^nmyt) 

= J2 PnroPm exp {-\a 2 (E n - E m )H) (108) 

n,m 

= J2 PnroPm exp {-\a 2 {E n - E m ft) + ^ P n f P n . 

Therefore, by inverting the transformation (p9|) , we obtain the solution of the Lindblad 
equation in the original Schrodinger picture as: 

Pt = E PnPoP m e- iiE "- E - )t -S' 72{En - Em)2t + Y, PnPoPn- (109) 

We recover the initial state po by setting t = in the right hand side of (|109|) . The 
off-diagonal terms are damped away exponentially at the rate la 2 V nm as t — > oo, where 
V nm = (E n — E m ) 2 is the square of the spread between relevant energy levels, and we 
are left with the Liiders state (^) for the terminal density matrix p^. 

It is interesting to observe that the Liiders state, obtained by the limit as t — > oo 
of the density matrix associated with the reduction process (11), coincides with the 
asymptotic time average of the density matrix in the case of a purely unitary evolution 
governed by the von Neumann equation 

^ = -i[H,p t ], (110) 

for which the solution is p t = e~ lHt f>oe lHt . More specifically, if we write 

{p)T±7 t £fH& (111) 

for the time average of p t up to time T, then we find that 

lim(p) T = VP n p P n , (112) 

T — >oo ' • 

n 

where P n is the projection operator onto the subspace TL n C H spanned by the states 
of energy E n . This result can be verified directly by use of the resolution of the identity 
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([TOD. The calculation is as follows: 

(p)t = ^ E / Pne- i6t poe™P m dt 

m,n u 

= i^P n p Pm/ e-^-^Mt (113) 



/ J PnPoPn + 7^ 2_^ PnPoPm ( 



sin(^ nm T) . cos(^ ram T) - 1 






where u; nm = E n — E m . Therefore, in the limit T — > oo the off-diagonal terms drop out, 
and we recover ( |112| ). For a closely related result see |24 . 



9. Change of Measure 

We return now to the stochastic differential equation (|TT|) governing quantum state 
reduction with a view to gaining further insights into the nature of the resulting 
dynamics. We shall demonstrate in this section how a 'change of measure' technique 
can be used to solve ([Ll]) and thus, in effect, to construct an explicit unravelling of the 
Lindblad equation (|98D . The general problem of formulating an appropriate unravelling 
of the Lindblad equation in a given physical context is a matter of considerable interest 



in a number of areas of modern physics p5| ^6], ^7], ^, [29], 

We begin with the following remark. Let jX t and a t be bounded JF r adapted operator- 
valued processes on (Q, JF, P) with the property that for all s, t G [0, oo) the random 
matrices fi S) ft t , o~ s and a t mutually commute. Then the stochastic differential equation 

d\iP t )=fi t \A)dt + cr t \il> t )dW t (114) 

has the unique solution 

|V t ) = exp (j (fi s - \a 2 s ) ds + J & s dW s j |Vo>- (H5) 

Here we allow for the possibility that the initial state \i/jo) may be random. A 
straightforward application of Ito's lemma shows that ( |115p leads back to ( p.l4| ). In 



the case of the reduction process (|Tl~D, which is evidently of the form ( |114| ), we can write 
ft = -ii?-yfff-ff t ) 2 , a t = \a ( H - H t ). (116) 



It follows therefore that 

\if>t) = exp (-iHt - \a 2 j (h- H^ ds 

+\o f (H-H^&WA^). (117) 

This is still an implicit solution for \ipt), because H s = {ip s \H\ip a ). Nevertheless, as a 
consequence of ( |117| ) we see that the evolution of the state vector according to flTT| ) can 
be expressed in the simple form 

IV*} = u t RM, (us) 
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where the operator-valued process U t is defined by 

U t 4 exp f-iHt) , (119) 



and the operator-valued process R t is defined by 

Rt ± exp (\a I (h- H^j dW s - \o 2 j [h - H^j * ds\ . (120) 

We note that Ut is unitary and that Ut and Rt commute. The square of i?j, which 
we denote by M t , is an operator- valued martingale. The fact that M t satisfies the 
martingale condition E s [M t ] = M s is evident from the expression 

M t = exp (a j (A- F s ) dW s - \a 2 j (h - H^ ds\ . (121) 

In particular, if P n is the projection operator onto the subspace Ti n C Ti spanned by 
the states of energy E n , then we find that 

M t = Y,PnM nt , (122) 

n 

where M nt is given by ([731). We note that for each value of n the process M nt is an 
exponential martingale. 

Now suppose that G is an observable that commutes with the Hamiltonian H. 
Then for its expectation in the state \ipt) we have 

G t = {A\G\A) 

= (^o\RtU}GUtRt\^o) (123) 

= (Vo|GM t |Vb), 
and therefore 

E s [G t ] =E s [(^ |GM t |Vo) 
= (^ \GE s [M t ]\^ ) 

= {MGM S \%) (124) 

= (ip \R s GR s \ip ) 

= {MR s utGU s R s \%) 

= (HG\lps): 

which shows that G t is a martingale. In this way we are able to verify directly that the 
dynamical law (]TT|) implies that the expectation value of any observable that commutes 
with the Hamiltonian is a weakly conserved quantity. 

To proceed further we note that it is a straightforward algebraic exercise to verify 
that M t can be expressed as the following quotient: 



exp (a C H(dW s + aH s ds) - \a 2 f* H 2 ds 

M t = f '-. (125) 

exp (a jl H s (dW s + aH s ds) - \o 2 /* H^ds 
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In particular, let us define the 'modified' Brownian motion process W^ by 

W* =W t + a I H s ds, (126) 

Jo 

so dWf = dWt + Hfdt. Then, because H is constant, we can write M t in the simple 
form 



M t = ij exp [aHW* - \a 2 H 2 i) , (127) 



where 



A* 4 exp (a ! H s dW* - \o 2 f H 2 ds\ . (128) 

The significance of the processes W£ and Aj will become apparent shortly. 

We have already verified that (pLlf) preserves the norm of |0o)- If we assume that 
("00 1 "00 ) = 1) then it follows from (|118|) that ( < o |^t|"0o) = 1 f° r ah £• Thus we deduce 



from (11251) and (|126|) that 



A t * = (0o I exp (aHW* - \a 2 H 2 t) |Vo)- (129) 



As a consequence we can write 

exp (aHW* - \a 2 H 2 t) 
Mt = y~„ ^ , (130) 

which has the important effect of localising the dependence of M t on H t in the modified 
Brownian motion W t *. The process H t in turn is given by (|T2"|), from which it follows 
that H t = (ip \HM t \ip ). Therefore, by use of ( |130|) we have 



(0o|#exp (aHW* - \a 2 H 2 t) |0 O ) 
Ht = ^- -^ , (131) 



V u 2 

which shows that H t can be expressed as a function of W t * and t. This is given explicitly 
by 

H T, n ^nE n e^{aE n W:-\a 2 E 2 n t) 
4 Z n Knexp{aE n W t *-^E*t) ' 

where as usual ir n denotes the probability that the eigenvalue attained is E n , given the 
initial state |0o). We also note that 

A* = J2 Kn exp (aE n W t * - \a 2 E 2 n t) , (133) 

n 

and that 

Now we proceed to examine the processes W£ and Aj* more closely. This is the point 
at which we introduce the highly useful concept of a change of probability measure. We 
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shall see in what follows that there exists a change of measure P — > Q such that for 
any given interval of time [0,T] the process W% for t G [0, T] is a Brownian motion 
with respect to the probability space (Q, Tt, Q) and the filtration Tt (0 < t < T). The 
implication of this is that with respect to the measure Q the basic processes H t , A£, and 
M t can be expressed in terms of ratios of sums of geometric Brownian motions, thus 
offering a significant element of analytic tractability. 

We begin with a few mathematical preliminaries concerning the change of measure 
technique. Given the probability space (fi, T, P), we recall that 1a denotes the indicator 
function of the event A G T. Thus for each uo G Vt we have Ia(^) = I if uj £ A and 
1,4.(0;) = if uj <£ A. It follows that 

Prob p [A]=E p [U], (135) 

where Prob and E p denote probability and expectation with respect to the measure P. 
Now let A be a positive random variable on the probability space (Jl,J-,¥). Then 
we can define a new probability measure Q on the underlying measurable space (Q, T) 
by the formula 

Pr°b«[.4] = ^f 1 - (136) 

Because A is positive, this relation is invertible and we have 

^-^ <-> 

where A* = 1/A. The two probability measures P and Q in this case are said to 
be equivalent in the sense that they agree on null sets, i.e., for all A G T we have 
Prob p [A] = if and only if Prob Q [A] = 0. 

In the case of a filtered probability space some important additional structure arises 
in this connection. Suppose the process A t is a positive martingale on (fl, J 7 , P) with 
respect to the filtration T, satisfying Ao = 1. For any fixed value of t the random 
variable A t can be used to define a measure Q t on (O, Tt) according to the procedure 
outlined in the previous paragraph. It follows then from (|136|) by virtue of the martingale 
property of A t that 

Proh Qt [A] =E v [A t l A ] (138) 

for all A G T t . We note that if s < t then Pioh Qs [A] = Proh Qt [A] for all A G T s . This 
is because 

Pioh Qt [A] =E p [A t l A ] 

= E p [E p [A t l A \T s ]] 

= E p [E p [A t |^]U] (139) 

= E P [A S U] 

= Pioh Qs [A\. 

Therefore, for any finite interval of time [0, T] the measure thus obtained on (Q, Tt) is 
independent of the specific choice of T. Thus we can drop the suffix on Q and speak of 
the change of measure P — > Q induced by the given density martingale A t . 
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The key result making use of this apparatus that we require in what follows is the 
theorem of Girsanov (see, e.g., [ftl|). Let [0, T] be a fixed interval of time, and W t a 
Brownian motion on (Q, Tt) with respect to the filtration Tt (0 < t < T) and the 
measure P. Suppose that the process \ t is jF r adapted and that 

A* = exp f- J \ s dW s - \ J \ 2 s ds\ . (140) 

is a martingale. Then Girsanov's theorem states that the modified process 

W *A Wt+ J XAs (m) 

Jo 
is a Brownian motion with respect to the equivalent measure Q induced by the density 
martingale A t . For A t to be a martingale it suffices that X s should satisfy the Novikov 
condition 

< oo. (142) 



E LJ 



exp(|^ A s 2 ds) 



In particular, if \ t is bounded, then A 4 is a martingale. 

If A* is a P-martingale then the associated process A* t = 1/Af given by 

A t * = exp (+ J \ s m: - \ J X 2 s ds) , (143) 

is a Q-martingale, and induces the inverse change of measure Q — > P. In particular, for 
any ^-measurable random variable X t we have the following formulae for the calculation 
of expectations: 

<[Xt] = iE? [A* t X t ] , (144) 

8 

and its reversal 



E?[X t ] = ^-Ef[A t XJ . (145) 

Returning to the matter at hand, we note that for quantum state reduction the 
process X t is given by o~H t , and the corresponding change of measure density martingale 
is given by 

A t = exp (-a I H s dW s - \a 2 I H 2 s ds\ (146) 

for the transformation P — ► Q. The process W* as defined by ( |126| ) is a Q-Brownian 
motion. The associated inverse transformation Q — * P is induced by the process A 4 
defined in (ggg). 

Now we are in a position to give a complete characterisation of the solution of 
the dynamical equation (|ll]) for the state reduction problem valid over any finite time 
interval [0, T] . The recipe is as follows. 

We start with the measure Q for which W* is a Brownian motion. Given W t * we 
then construct the process H t by use of formula (|132|) , and the process A*[ by use of 
formula Ql33|) , and the process M t by use of formula (|134j) . Thus we see that the wave 
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function \ip t ) along with all the related processes R t , M t , A t , and H t can be explicitly 
constructed as functions of W£ and t. The physical measure P constructed by use of 
A* is then used for the calculation of ensemble averages. In particular, letting E denote 
the expectation with respect to the physical measure P, it follows from ( |134| ) that 

E [X t ] = E Q [A*X t ] (147) 

for any jF t -measurable random variable X t . 

For example, suppose G is an observable that does not necessarily commute with 
the Hamiltonian H, and we wish to calculate the ensemble average of the expectation 
value (ip t \G\ip t ). Then by use of (|147|) we have 



E 



<^|G|Vt> 



E^ 



E^ 



E c 



KWt\G\rl>t) 

A* t {MUlRtGR t Ut\^) 

(V'ole 



148) 



iHt+\uHW*-\a 2 H 2 tQ -Mt+\aHW*- 



\a 2 H 2 t 



l^o) 



V" G pKEm-E^t+laiEm+E^Wf-^El+E^t 



where the matrix elements G mn are given by 

G mn =(iP \P m GPM- (149) 

Since W^ is normally distributed with mean zero and variance t with respect to the 
Q-measure, the expectation in (|148|) can be readily computed. By use of the simple 
relation 



which holds for any constant a, we see that 

E Q [ e la(E m +E n )W*-\(El+El)a 2 i 



\a 2 (E m -E n ) 2 t 



As a consequence we deduce that 



where 



E 



G, 



(&|G|&>1 =Y, r. 



j(E m —E n )t—-za (E m —E„) t 



(150) 
(151) 
(152) 



E 



(153) 



(iPo\P m GP n \ip ) 

= Tip (P m GP n ) 

= TrG(P n p P m ) 

and po is the density matrix corresponding to the random initial state. This result is 
consistent with our earlier expression fll09| ) for the solution of the Lindblad equation, 
and illustrates the fact that the change of measure technique is indeed highly effective 
as a calculational tool for quantum state reduction models. 
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